Mutual Impedance of Grounded Wires Lying on the 
Surface of the Earth * 

By RONALD M. FOSTER 

This paper presents a formula for the mutual impedance between two 
insulated wires of negligible diameter lying on the surface of the earth and 
grounded at their end-points. The formula holds for frequencies which are 
not too high to allow all displacement currents to be neglected. For any 
two elements dS, ds of the two wires the mutual impedance is obtained 
from their direct-current mutual impedance by introducing the complex 
factor 2 (yr)--[_l - (1 + yr)«— >*] in the reactance term, y being the propaga- 
tion constant in the earth, and r the distance between the elements dS 
and ds. 

THE mutual impedance of grounded circuits may be derived from 
certain results obtained by A. Sommerfeld, 1 who has developed 
formula? for the electric and magnetic fields in the earth and in the air 
due to horizontal and vertical electric and magnetic antennae situated 
at the surface of the earth. For our present problem we use his formu- 
lae for the electric field in the earth due to a horizontal electric doublet, 
since this doublet may be regarded as a short element dS of a wire of 
negligible diameter carrying a finite current. At the end of this 
present paper we shall show how the same formula for the mutual 
impedance may be obtained directly from first principles. 

Sommerfeld uses rectangular coordinates (x, y, z) and the corre- 
sponding cylindrical coordinates (r, <f>, z), the surface of the earth, 
assumed flat, being the xy plane, and the z axis extending upward into 
the air. The doublet is at the origin, and its axis along the x axis. 
Then the components of the Hertzian vector 2 in the earth (z < 0) from 
which the electric field is determined are 3 

(1) v.- c% £*&***+, 

(2) n„ = o, 

* Presented by title at the Eugene, Oregon meeting of the American Mathemat- 
ical Society, June 20, 1930, as "Mutual Impedances of Grounded Circuits." 

l A Sommerfeld, " Ober die Ausbreitung der Wellen in der drahtlosen Tele- 
graphic" Annalen der Physik, (4), 81, 1135-1153 (December 1926). This paper is a 
summary and an extension of earlier work by Sommerfeld and von Hoerschelmann, 
references to which will be found in the paper. ... . 

*H. Abraham and A. Foppl, "Theorie der Elektrizitiit, ' 5th ed., Leipzig and 
Berlin,' 1918; Vol. I, § 79, page 331. 

3 A. Sommerfeld, loc. cit., pages 1145 and 1146, introducing the constant lactor 
defined on page 1152. 
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IT, = CW - *o 2 ) j£ cos f" ^p- e'^fdp, 



where the time factor e~ iul is omitted throughout. Jo is the Bessel 
function of order zero, and the constants k and ko are the propagation 
constants in the earth and in the air for plane waves varying with the 
time as c~ iut . Their values in Heaviside units are given by Sommerfeld 
as 

(4) k- = -j (ew 2 + iau), k 2 = -s e w 2 , 

where e and e are the dielectric constants of the earth and of the air, 
respectively, a is the conductivity of the earth, assumed uniform, and 
c is the velocity of light. In both media the permeability is taken as 
unity. Also 



(5) N = 1H? - k * + feoWp 2 - #, 



(6) N' = Vp 2 - k - + Vp 2 - k\ 

and C is a constant measuring the electric moment of the doublet. 

We now replace the doublet by a short element of wire dS carrying a 
current Ie i(ot , and at the same time we assume that e and e are both 
negligible, so that all displacement currents are neglected. This is a 
simplification which is ordinarily made as a first approximation at 
power frequencies for the shorter transmission lines. Then, introduc- 
ing c.g.s. electromagnetic units, in which the conductivity of the 
earth is X, and noting that we have changed the sign of co, formula? 
(4)-(6) become 

(7) k- = — 747rXco = — 7 2 , 

(8) £o 2 = 0, 

(9) N= - t 2 P, 

(10) N' = P + Vp 2 + i\ 
and the constant C is such that 

(11) -rr- = - X current X effective length of doublet 

H (S 

= IdS 
2ttX' 
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Substituting from (7)— (11) in (l)-(3) we have, therefore, 

(12) n,-gfr J t p \ ***** 

2 ^ Jo P + V7T7 

IdS ( d 2 P d 3 Q d 3 Q 



27r\7 2 \ds i dx*dz dy 2 d 

(13) n„ = o. 

7</5 / &P d a Q 

2 
where 



2ir\y- \ dxdz dxdz 2 



(15) P= r/o(rp)6^"^-=' 

Jo Vp- 



Vp' J + 7 2 



= — e -y« 



and 
(16) 



Jo Vp 2 + 7" 

= j fl [| 7 (^ + i)Wh(« - »)1 



with i? 2 = r 2 + s 2 . 

The integral P is well known, 4 while () is evaluated by a suitable 
transformation of a Fourier integral. 5 h(z) = Jo(iz) and K (z) 
= $ir*2J (l) (*s) are the Bessel functions of the first and second kinds 
for imaginary arguments as denned by G. N. Watson. 6 In reducing 
IIx to this form we use the differential equation 7 for J to obtain the 
relation 



(£+^) j ^ + " u ^' - 



The components of the electric force in the earth are obtained from 
II by the formula 



(17) E = grad div II - ? 2 II, 

«See 
page 72 



4 See e.g. H. Bateman, "Electrical and Optical Wave-Motion," Cambridge, 1915, 
; or G. N. Watson, "Theory of Bessel Functions," Cambridge, 1922, page 416, 



formula (2) of § 13.47, with /* = and v = -J. 

6 G. A. Campbell, "The Practical Application of the Fourier Integral, Bell 
Svstem Technical Journal, 7, 639-707; using pair 936 of Table I, with or = 5. substitut- 
ing x 2 for (g 2 - 4) in the integral of G, and generalizing the resulting integral to in- 
clude complex quantities. 

■ G. N. Watson, op. tit., pages 77, 78. 

7 G. N. Watson, op. cit., page 19, formula (1) of § 2.13. 



MUTUAL IMPEDANCE OF GROUNDED WIRES 411 

and we thus obtain E x , E y , E : in the compact form 

nn (f f f \ - IdS ( - i!&. - d2p g!g d * p \ 

1 j ■ " "' z) 2ttXV 3/da ai 5 "' d.ro^s' any- 

where P and (J are given by (15) and (16). In deriving this form we 
use the fact that Q satisfies the wave equation 

^ + ^ + 0.^=0. 

dx 2 By- dz 1 

At the surface of the earth (z = 0) the electric force takes the simple 
form 



<■« ^-m-uv +i 



+ yr _ yr a ? / 1 
— c yr 



r 1 dxdy \ r 

where we have used the expressions for the derivatives 8 of the Bessel 
functions, I '(z) = h{z), K '(z) = - Ki(z), and also the identity 9 
h(z)Kt(z) + h{z)K Q (z) = 1/e. 

The mutual impedance dZ n between two infinitesimal elements dS 
and ds is now written down as the ratio of the resulting electric force in 
one element to the current in the other, with sign reversed: 

dSds\ cV- /1\ \ + yr 

(20) rfZ I2 = — ^ose-^j-cose-^-,- 

■ d 2 /I 

dxdy \ r 



— sin e 

dSds f 3 sin $ sin — cos e cos e 

2ttX 

2 



[3 sin 4> sin (/> — cos e cos e . . . I 

[ ^ ^-(! + yr)eryrl 

Sds I d 2 / \\ , cos « r , ,. , . _1 



where * and are the angles which the elements dS and ds make with 
r, and e = <f> — is the angle they make with each other. 

Integration over the two wires 5 and 5 gives a general formula for 
the mutual impedance of grounded wires lying on the surface of the 
earth : 



(21) 






1 d- I 1 

2tt\ ' dSds 



iu C ~^r{ Wr [1 - (1 + Tr)c " 7r] } ] dSds - 



8 G. N. Watson, op. cil., page 79, formula (7) of § 3.71. 

G. N. Watson, op. cit., page 80, formula (20) of § 3.71, with v = 0. 
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The factor 
(22) 



(7') ! 



[1 - (1 + yr)(T»] 



approaches unity as w approaches zero, and Zu then agrees with the 
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Fig. 1— Real and imaginary parts of the complex factor, 
plotted as functions of r' "\yr\ = (iir\w) ll2 r. 
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direct-current mutual impedance as given by G. A. Campbell. 10 
Introducing this factor, which is a function of yr only, into the re- 
actance term for the direct-current mutual impedance between two 
elements dS and ds gives the general expression for their mutual im- 
pedance corresponding to the propagation constant 7. It is interesting 
also to determine, for any given value of 7, the variation of the factor 
(22) for increasing values of r. This is shown very clearly in Fig. 1, 
where the real and imaginary parts of (22) are plotted for increasing 
values of r' = \yr\ = (47r\oj) l/2 r. The real part, we note, decreases 
rapidly from the initial value unity as r' increases, while the imaginary 
part is always negative, decreasing from zero to a minimum value 
(approximately — 0.3 for r' = 1.5) and then increasing towards zero, 
although it does not approach zero so rapidly as the real part does. 

The first three terms in the expansion of Z 12 for low frequencies 
are given by 

+ (1 - *)H87rW)' '- A B ab cos 6 + ■■•, 

where N Ss is the mutual Neumann integral between the two wires 5 
and s of arbitrary form but with end-points A, B and a, b respectively; 
6 is the angle between the straight lines AB and ab. The first two 
terms in this expansion are precisely the direct-current mutual im- 
pedance as given by G. A. Campbell. 

The first term in the expansion of Z i2 for a long straight wire S and 
any wire 5 located near the midpoint of 5 is 



(24) 



j -7-9 z~Ki(yx) coseds 

J L TA.V- 7TA.V 



x being the positive distance from ds to S, and e the angle between ds 
and S. Ki(z) = — %TlIi (l) (iz) is the Bessel function of the second 
kind for imaginary argument as defined by G. N. Watson. 11 In ob- 
taining (24) from (21) we use the derivative with respect to x of the 
integral 



f 



e~ yr 

— <lz = K (yx), 



which is a special case of the integral used above in evaluating Q, with 
x assumed positive. 

10 G. A. Campbell, "Mutual Impedances of Grounded Circuits," Bell System 
Technical Journal, 2, 1-30 (October 1923). 

11 G. N. Watson, op. cit., page 78. 
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The expression in square brackets in (24) is the mutual impedance 
gradient parallel to an infinite wire at a positive distance x from the 
wire. It agrees with the results published independently by F. 
Pollaczek, 12 J. R. Carson, 13 and G. Haberland, 14 and has been employed 
by us to obtain numerical results since 1917. Pollaczek has also in- 
vestigated the case of two gounded circuits of finite length. 15 

The mutual impedance dZ n between a short grounded circuit dS 
and a counterclockwise small loop of area da, on the surface of the 
earth, is given by the formula 

(25) dZ 12 = *S5 • S -^ [3 - (3 + Syr + y^)^J 

where <p is the angle which dS makes with r, the line from da to dS. 
This may be obtained from Sommerfeld's formulae for the horizontal 
electric force due to a vertical magnetic antenna, or it may be obtained 
by an application of Stokes's theorem to formula (20) above. 

By a further application of Stokes's theorem we may obtain the 
mutual impedance between two counterclockwise small loops dA and 
da, namely, 

(26) dZ 12 = ~ • \ [(9 + 9yr + 4 7 2 r 2 + y^)e^ - 9]. 

2. 7TA T 

This result might also be derived from Sommerfeld's formula for the 
vertical magnetic force due to a vertical magnetic antenna. 

We shall now indicate briefly how the same value of E as given in 
(18) above may be obtained directly, though more laboriously, from 
first principles. In this method we start from the fundamental 
solution 18 

(27) u = e lx+mv+n *e iat 
of the wave equation 

12 F. Pollaczek, " Uber das Feld einer unendlich langen wechselstromdurchflossenen 
Einfachleitung," Elektnsche Nachrichten-lechnik, 3, 339-359 (September 1926). 

13 J. R. Carson, "Wave Propagation in Overhead Wires with Ground Return,' 
Bell System Technical Journal, 5, 539-554 (October 1926) ff 

14 G. Haberland, "Theorie der Leitung von Wechselstrom durch die krde, 
Zeitschrift fur angewandte Mathematik und Mechanik, 6, 366-379 (October 1926). 

16 F Pollaczek " Gegenseitige Induktion zwischen Wechselstromfreileitungen von 
endlich'er Lange," Annalen der Physik, (4), 87, 965-999 (December 1928). His as- 
sumptions regarding conditions at the ground connections seem to depart considerably 
from the conditions assumed in the present paper, and moreover his results are not 
expressed in convenient form for direct comparison with the formula given above 

16 H. Bateman, op. cU.,\ 4, pages 6, 7; §11, page 26. 
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which is satisfied by the electric force in the earth; 7 = (i47rXw) l/2 is 
the propagation constant for plane waves which vary with the time as 
g** 1 . The parameters /, m, n satisfy the relation 

(29) I 2 + m 2 + n 2 - y 2 = 0. 

In the air, the same equations hold, but with the propagation constant 
7 equal to zero, and we note that the solution in the air must be chosen 
to vanish at an infinite height, while in the earth the solution must 
vanish at an infinite depth. 

For convenience in this method we start with a short straight wire of 
length 2a lying along the x axis, later allowing a to approach zero. 
Thus we suppose that the current Ie iut enters the earth at the point 
(a, 0, 0) and leaves it at the point (-a, 0, 0). The factor e ut will be 
omitted, however, throughout the following work. The current flow 
in this system is symmetrical with respect to the vertical plane through 
the wire, the xz plane, and is also symmetrical, but with sign reversed, 
with respect to the vertical plane normal to the wire at its midpoint, 
the yz plane. Then if we replace the three parameters /, m, n of (27) 
by two independent parameters ju, v, such that 

(30) / = ± in, m = ± iv, n = ± Vm" + v* + T, 

formula (29) is identically satisfied, and we can then replace the four 
solutions £**'***«' by their corresponding expressions in terms of sines 
and cosines, namely, 

sin xn sin yv, sin xp cos yv, cos x/t sin yv, cos x/i cos yv. 

The above considerations of symmetry will eliminate, for each com- 
ponent of the electric force, all but one of these forms. With the re- 
maining solution as a basis we build up, by means of the Fourier in- 
tegral, a general expression for any possible steady harmonic oscilla- 
tion. Hence we may write down the general solutions for the total 
electric force in the earth (2 < 0), as follows. 

I F x (n, ^e^+'H-^cos xn cos yvdpdv, 
Jo 

/"<"> /•* 

(32) E u = I F y (p, v)e z ^ 2 +"+-> 1 sin x/x sin yvd/xdv, 

Jo t/0 

J '•CO /'•CO 
I F,{p, ^MH-'+T'gjn Xfl cos y v did Vl 
»/0 
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where the positive sign is chosen in the exponential term containing z 
since the solution must vanish at an infinite depth, z being negative in 
the earth ; and that value of the radical is taken which has a positive 
real part. F x , F v , F z are arbitrary functions of their arguments, to be 
determined by the physical conditions of the problem. 

In the air (0 < z) we may formulate the corresponding solutions for 
the total electric force as 

(34) E x = I P x (n,v)e- zy[] *+ ii cos xn cos yvdndv, 

Jo Jo 

/-•CO /»co ^_^_ 

(35) E v = I I P v (fi, v)e-™"+" sin xii sin yvdftdv, 

Jo Jo 

/-»X /»x 

(36) E z = I I P t (n, v)e~" /i5 + ;i sin Xfi cos yvdfidv, 

Jo Jo 

where the propagation constant is zero in the air; the negative sign is 
chosen in the exponential term containing z since the solution must 
vanish at an infinite height, z being positive in the air; and P x , P u , P-. 
are arbitrary functions of their arguments. 

To determine these six arbitrary functions we need six independent 
relations among them. Two of these relations are obtained by 
utilizing the fact that the divergence of the electric force either in the 
earth or in the air is equal to zero, that is, 

BE, dEy dE, = 
dx ^ dy "•" dz 

By means of this we obtain from (31)-(33), 

(37) - nF x + vF v + Vm 2 + ? + V*F* = 0. 
and from (34)-(36), 



(38) - iiP x + vP y - a/m' + *P, = 0. 

Since the horizontal components of the electric force are continuous 
at the surface of the earth (z = 0) we see that we must also have, from 
(31) and (34), 

(39) F x = P x , 
and from (32) and (35), 

(40) F v = P v . 
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We may obtain a fifth relation from the fact that the current / 
flows through the earth from one grounding point to the other. To 
utilize this fact let us compute the total current flowing out through 
five faces of a rectangular prism in the earth, the sixth face being a 
rectangle in the surface of the earth surrounding the grounding point 
(a, 0, 0), the prism extending from x = a — £ to x = a 4- £, from 
y = — i\ to y = 7], and from z = — f to s = 0. The components of 
the electric force being given by (31)-(33), and X being the conductivity 
of the earth, we obtain for this current the expression 

(41) - 4X f r^suiaMsinfMsin,, 

Jo Jo V" 

after simplifying by means of the divergence condition (37). This 
current flowing out through the prism is / if the face in the surface of 
the earth includes only the one grounding point (a, 0, 0), but is zero 
if it includes both grounding points; that is, the above integral (41) 
equals / if £ < 2a, but equals zero if 2a < £, for any positive value of ??. 
It is readily verified that the Fourier integral 



(42) — f" f* sin* an sin 

^Jo X M" 



£ H sin tw 
ajxai 



has the desired properties. Accordingly, we must have 

21 

(43) F z = — ^sina^. 

7T-A 

To obtain the one additional relation which is needed, we make use 
of the fact that the current / flows through the straight wire from one 
grounding point to the other. Let us integrate the magnetic force 
around a rectangle in a plane perpendicular to the wire, that is, 
perpendicular to the x axis, the rectangle extending from y = — v 
to y — 77 and from z = — J* to z = f , the path of integration being 
taken in the clockwise direction looking along the positive direction of 
the x axis, and then equate this integral to 4tt times the total current 
threading the rectangle. The components of the magnetic force 
which we need, H y and H z , are found from the fact that curl E= -iuH, 
that is, 

(44) iuHv = fl- d lk t 

dx dz 

(45) faff. = *&_*& 

dy dx 
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where the E's are given by (31)-(33) for z < and by (34)-(36) for 
< z. We now subtract from this integral 4tt times the current in 
the earth which threads the rectangle, this quantity being found by the 
appropriate integration of E x , as given by (31), over that portion of the 
area of the rectangle which lies below the surface of the earth. As a 
final result we obtain the expression 

(46 ) 3_ f" f "i (- V M 2 + ? + y°-F x + nF, - Vm 2 + *P* - nP t ) 

w Jo Jo 

X cos xy sin rjvdydv, 

after simplifying by means of the divergence conditions (37) and (38). 
The net current threading the rectangle, after subtracting the current 
in the earth, is / if the rectangle is situated between the two grounding 
points, but is zero if it is outside them ; that is, the above integral (46) 
equals M if \x\ < a, but equals zero if a < \x\ , for any positive value 
of t]. It is readily verified that the Fourier integral 

167 ("* f"° sin ay cos xy sin yv , , 

~Jo Jo »* 

has the desired properties. Accordingly we must have 



(47) - Vm 2 + ? + -fF* + m/ 7 . - Vm 2 + *P* - nP. 

Si<al sin ay 



T H 



We can now solve equations (37)-(40), (43), and (47) for the six 
arbitrary functions, obtaining 

(48) F,= P, = ^r-=^-3^ + ^iX ! ]sina„ 
tt-X L m Vm 2 + v l y J 



27 v . 

(49) F v = P v = -5r • , sin a/t, 

21 . 
(43) P* = ?r sin ay, 



-X 



_ 2/ _ Vj n 2 + y 2 + t 



(50) P, = ^V • AM - 1 " " PT sin a "- 



Substituting these values in equations (31)-(33) and letting a 
approach zero such that la = dS, we find, for the electric force in the 
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earth, 

X cos xn cos yv dfidv, 

(52) E y = ^ r f X ; P " e*^^ 5 sin * M sin yp d-ndv, 

™ Jo Jo \M a + "* 

(53) E z = - -r- f f /xe iV ' l2+, ' 2+72 sin .r/x cos yv dpdv. 

™ J J 

These are precisely the values found by the former method, for the 
integrals P and Q may be expressed as double integrals by substituting 
for J (rp) the integral expression given by the formula 17 

(54) 7 (Wm 2 + "*) = - I cos(r/i cos 0) cos (r* sin 6)dd, 

* Jo 

and introducing rectangular coordinates in place of r, 0. These inte- 
grals may, therefore, be written in the equivalent forms, 



(55) P = — I I ms x [teas -vv dfidv, 

^Jo Jo Vm 2 + "' 2 + 7 2 



r VM2+f'+7 8 



9 /-»oo /'•oo ^V^+^+T* 

(56) = - I I , ,- , n cosxn cos yvdfidv, 

* Jo Jo V/x' J + "W + "- + 7" 

and comparison with (51)-(53) again leads to the values 

_ IdS I _ d 3 (? _ d^P gg d'-P \ 

(18) (£ x , £„, £.0 - y^ ^ dy2da as , . aA . a>)a2 ' a.vas/' 

where P and Q are evaluated in (15) and (16). Thus the mutual im- 
pedance formula presented in this paper may be derived directly from 
first principles, without reference to the work of Sommerfeld. 

I am greatly indebted to my colleague, Dr. Marion C. Gray, for 
putting into its present form the derivation of my formula from 
Sommerfeld's results. 

17 G. N. Watson, op. cil., page 21, formula (1) of § 2.21. 



